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The total spacetime manifold for a Schwarzschild black hole (BH) is described by the Kruskal 
coordinates u = u{r,t) and v = v{r,t), where r and t are the conventional Schwarzschild radial 
and time coordinates respectively. The relationship between r and t for a test particle moving on a 
radial or non-radial geodesic is well known. Similarly, the expression for the vacuum Schwarzschild 
derivative for a geodesic, in terms of the constants of motion, is well known. However, the same is 
not true for the Kruskal coordinates; and, we derive here the expression for the Kruskal derivative 
for a radial geodesic in terms of the constants of motion. In particular, it is seen that the value 
of I du/dv I (= 1) is regular on the Event Horizon of the Black Hole. The regular nature of the 
Kruskal derivative is in sharp contrast with the Schwarzschild derivative, | dt/dr \= oo, at the Event 
Horizon. We also explicitly obtain the value of the Kruskal coordinates on the Event Horizon as a 
function of the constant of motion for a test particle on a radial geodesic. 
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I. INTRODUCTION 



It is known for more than 80 years that the region exterior to a point mass or the event horizon {r > Vg — 2m) of 
a Schwarzschild Black Hole (BH) can be described by the vacuum Schwarzschild metric 

ds'^ = gttdt^ + grrdr"^ + gBedB"^ + g^4,d(f' (1) 

where gu = (1 — 2m/r), grr = —(1 — 2rn/r)~^, gee = — r^, and g^^ = — sin^ 0. Here, we are working with a 
spacetime signature of +1, -1, -1, -1 and r has a distinct physical significance as the invariant circumference radius. 
The coordinate time t too has a physical significance as the proper time of a distant inertial observer 5*00- At r = 2m, 
grr blows up and as r < 2m, the gu and grr suddenly exchange their signatures though the signatures of ggg and 
remain unchanged. The detail dynamics of a "test particle" in the vacuum external spacetime is well known for 
a very long time and discussion on it is contained in practically every text book or monograpgh on classical General 
Theory of Relativity (GTR). One of the key aspects for studying the kinematics of a test particle is the knowledge 
about the relevant derivative of the spatial coordinate with the temporal one. For instance for any geodesic having 
angular momentum or not, one knows the details about the behaviour of the Schwarzschild derivative dr/dt or dt/dr. 
And the fact that dt/dr = —oo blows up at the Event Horizon restricts the utility of the Schwarzschild dynamics 
below the EH. This tantamounts to the well known fact that the vacuum Schwarzschild metric fails to describe the 
spacetime inside r < 2m. 

On the other hand, we learnt in 1960 that both the exterior and the interior regions of a BH may be described by 
a one-piece coordinate system suggested by Kruskal and Szekeres |^,^. Though, in the intervening 39 years hundreds 
of articles have been written on Kruskal coordinates, and, most of the treatises on GTR too regularrly deliberate 
upon the original work of Kruskal and Szekers, the fact remains that sufficient effort has not been made to study 
the kinematics of a test particle in terms of the Kruskal coordinates so that one could have a better insight and 
appreciation of the kinematics inside the EH. As a first step towards this direction, in this paper, we would derive 
expressions for the Kruskal derivative du/dv for a radial geodesic. For the sake of completeness, we shall start from 
the usual description about the Kruskal coordinates and first derive the exact expression for the value of the Kruskal 
coordinates on the EH [uh and vh) in terms of r and t. We shall show here that uh and vh are always non-zero 
and finite in general. More importantly, we shall explicitly show that unlike the Schwarzschild derivative, the Kruskal 
derivates are regular on the EH in accordance with the singularity free nature of the Kruskal coordinates. 



II. KRUSKAL COORDINATES 

For the region exterior to the EH (Sectors I & HI), the Kruskal coordinates are defined as follows: 

u = /i(r)cosh- — ; v — fi(r) sinh - — ; r > 2m (2) 
4m 4m 

where 

A(r)=±(-^-l)'^%'-/4'" (3) 
v 2m / 



Here the plus sign corresponds to "our universe" while the negative sign corresponds to the "other universe" . The 
"other universe" is a legitimate mathematical solution of the Schwarzschild problem (irrespective of its observational 
reality), and is a time reversed mirror image of "our universe". 

And for the region interior to the horizon (Sectors II & IV), we have 

u = f2 (r) sinh - — ; v — f2 (r) cosh - — ; r < 2m (4) 
4m 4m 



where 



/2(r)=±(l--^V^%'-/4™ (5) 



2m/ 

In terms of u and v, the metric for the entire spacetime is 

r 

The metric coefficients are regular everywhere except at the intrinsic singularity r = 0, as is expected. Since afterall 
the Kruskal coordinates are defined using r and i, for a proper understanding of the Kruskal dynamics, it is necessary 
to recall the inter-relationship between the Schwarzschild coordinates for a geodesic. 
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A. Inter Relation Between Schwarzschild Coordinates 



For a test particle on a radial geodesic, the angular momentum is zero, and there is only one conserved quantity, 
the energy of the particle (per unit rest mass), E, as measured by a distant inertial observer: 



E = — 1 - 2m r) 
as 



(7) 



where s is the proper time. For a massless particle like a photon, we have E — oo, otherwise E is finite. For a radial 
geodesic, the motion of the particle is determined by (see Chandrasekhar, pp. 98) 1^ 



dr 
ds 



and 



so that 



— = -^E^-{l-2m/r) 

dt _ E 
ds 1 — 2m /r 

E{1 ~ 2m/r)-i 
v/-B2 _ (1 _ 2m/r) 



dt 
dr 



(8) 



(9) 



(10) 



Clearly as r ^ 2m, dt/ dr — > — oo. On the other hand, we do not expect such irregular behaviour for the Kruskal 
derivative. 

Here note that if the particle is released from rest [dr/ds = 0) at r = at t = 0, from Eq. (8), it is seen that j|] 



E^ ^{l- 2m/ri) 



or, 



r,/2m = (1 - £;2)-i 
It is convenient to introduce a (cylic) parameter 77 through 

„2 



r = ^{l + cosT]) = y— ^cos^(7?/2) = r,cos^(?7/2) 



Obviously, 77 = when r = ri and at the EH, we have 

V = Vh ~ 2 arcsini?; r ~ 2m 
Now after some manipulation, Chandrasekhar arrived at the following Eq. involving t and ?/ : 



'±^E 
dr/ 



2m 



1/2 



cos4(?7/2) 



^(r//2) -cos2(r/H/2) 



(11) 
(12) 

(13) 
(14) 

(15) 



This Eq. can be integrated to find the exact relation between t and r for a radial geodesic (actually, even for non-radial 
geodesic this Eq. would hold good): 



t 
2m 



E ii^)"^ + sin 77) + (1 - E-^U + In 



\ 2m 7 



tan(7]_H-/2) + tan(?7/2) 
tan(77H/2) - tan(7;/2) 



(16) 



The above Eq. may also be written without introducing t/h and E explicitly: (see pp. 824 of ref.[l] or pp. 343 of 
ref.[2]): 



t 1 I (r,/2m- 1)1/2 + (r,/r- 1)1/2 1/2 / x 

^=^'' " (r./2m-l)V2, 1)1/2 I + - 1 j h + [-^) + «-^) 



(17) 



We find from Eqs. (16-17) that, as r 2m from Sector I, the logarithmic term blows up and t — > 00, which is a well 
known result. Further Kruskal coordinates envisage that approach to the EH from the Sectors HI & IV corresponds 
to t = —00. 
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B. Kruskal Coordinates on the Event Horizon 



In Sectors I & III, Kruskal coordinates obey the relation 



u , t 
— = cotn - — 
V 4m 



And since r 2m corresponds to t ^ ±00, at the EH, we have 

r = 2m 



Vh 



On the other hand, in Sectors II and IV, we see 



u , t 
— = tanh - — 

V 4m 



(18) 



(19) 



(20) 



and as r ^ 2m, t ±00, we are led to the same Eq. (19). In the same limit, r 2m and t ±00, we find that 



Uh = Vh^ fi exp — 



(21) 



It might appear that since /i(2m) = /2(2m) = on the EH, we would have uh = ±vh = 0. But this is incorrect 
because the temporal part of u and v tend to blow up much more rapidly on the EH. And one has to carefully obtain 
the actual values of uh and vh by working out appropriate limits 
To do so we introduce a new variable 



z = ri /2m — 1 = 



and, let, in the vicinity of the EH, 



so that 



r/2TO = 1 + e; e ^ 



Then, in the vicinity of the EH, by retaining terms first order in e, we can rewrite Eq.(17) as 

. ^1/2 I ^1/2 _ . sl/2r .r\ 



z-,^ - z^i- (1 - 

V Amz) 

As e ^ 0, the logarithmic term in the above expression becomes 

1 + 1 



A{t) = In I 



4mz 



1-1 + 



In- 



8mz 



4mz 



Then, using Eqs. (24) and (26), we find 

/i'exp(A) 



8emz 



4e(l - 2m/ri) = 4eE'^ 



Now considering the other terms in the expression for t/2m in Eq.(25), we find that, in this limit, 

1/2 



2m 



u% = v% = 4e(l — 2m/ri) exp 
In terms of E, we have 

u]j=v]j= AeE'^ exp E 



1 



riH + 



r]H + sm rjH 



2Vl - E^ 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 
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One would have uh = vh = if = or, if the test particle is injected from, rest right at the EH. Clearly, this is 
unphysical, and thus we see that uh and vh are non-zero. Further, for a finite value of ri/2m oi for E < 1, they are 
finite too. The finiteness of u and v at the EH is physically appealing because u and v are expected to be completely 
regular at the EH. However for ri/2m = oo or _B = 1, we find = vjj = oo. 

On the other hand, since r = rj at f = 0, by using the definition of u and v, we find that the initial values of 

u^=u^ = {n/2m-l) = ^^ (30) 

and 

=vf = (31) 



III. KRUSKAL DERIVATIVE: A DIRECT APPROACH 



Having shown that uh and vh are, in general, non-zero, we are now in a position to evaluate the Kruskal derivative, 
the key ingredient for studying the Kruskal dynamics for a radial geodesic. We first confine ourselves to Sector I. By 
differentiating /i(r) (Eq.[3]) with r we obtain 

|^ = ^^(r/2m-l)-V^ (32) 
dr 2m 4m ^ ' ^ ^ ' 

Then by directly differentiating Eq.(2) by r, we find that irrespective of the sign of dfi/dr, we will have 

du dfi ^ t , h ■ t dt 

-r = -T cosh 1- — smh — — (33) 

dr dr Am 4m 4m dr 

Interestingly, in all the sectors, we obtain the same functional form of du/dr. Using Eqs.(2) and (4) in the foregoing 
Eq., we see that 



du u , , N 1 V dt 

3- = 1 - 2m/r)-i + 1— -r 

dr 4m 4m dr 

On the other hand by differentiating Eqs. (4) and (5), we find that 



-rr pr/im 



df2 



dr 2m 4m 



(r/2m- 1)- 



and 



du df2 . , t f2 t dt 

-r = ~r Sinn \- — cosh - — — 

dr dr 4m 4m 4m dr 



(34) 



(35) 



(36) 



And by using Eqs. (4) and (35) into the foregoing Eq., we obtain the same expression (34) for du/dr in Sectors H & 
IV. Further, using Eq.(lO) in (34), we obtain the ultimate expression for 



du 2m/r)-i 
dr 4m 



vE 



^JE^ - 1 + 2m/r 



(37) 



valid in all the sectors. Similarly, we obtain the ultimate functional form of dv/dr which is valid for all the sectors: 



dv 
dr 



(1 - 2iii/ry 
4m 



iiE 



- 1 + 2m/r 



(38) 



And, the general value of du/dv in any Sector is obtained by dividing Eq.(37) with (38): 



du 
dv 



^ E'^-l+2m/r 



uE 



(39) 
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A. Kruskal Derivative at the Event Horizon 



Since u and v are expected to be differentiable smooth continuous (singularity free) functions everywhere except at 
r = 0, and also since the "other universe" is a mirror image of "our universe", wc expect that the value of du/dv for 
any given r must be the same, except for a probable difference in the signature, in both the universes. The meaningful 
way to find the value of du/dr at the EH will be to concentrate on the Sectors II & IV for which uh = —vh 

du u-v 2uh ^ „ 

> = = —1; r = 2m (40) 

dv V — u 2vh 

The Eq. (39) for the Kruskal redivative, however, tends to yield a "0/0" form at r = 2m for Sectors I & III having 

Uh = Vh- But as mentioned above, we expect this 0/0 form to acquire the value du/dv = +1 because these Sectors are 
the mirror images of Sectors II & IV. Otherwise the whole idea of having an extended time symmetric Schwarzschild 
manifold would be inconsistent. Thus, in general, we must have 

du/dv = ±1; r — 2m (41) 

The fact that we must have du/dv = +1 for the Sectors I & III can be reconfirmed in the limiting case of u\ = vjj = oo 
for E = 1 or uh = vh = for the (unphysical case) E = directly by using L' Hospital's theorem. 
Note that, by this rule, we can write, 

u du/dr 
lim - = Hm — '— (42) 

u^O(oo),-y^O(oo) V u->0(oo),v->0(cx3) dv / dr 

In any case, from Eq. (19), we already know that u/v = ±1 at r = 2m. Then we can rearrange the foregoing Eq. 

as 

lim ^ = ±1 (43) 

■u^o(oo),t)^o(oo) dv/dr 



or. 



lim ^=±1 (44) 

r— >2m dv 



IV. A DIFFERENT ROUTE 

It may be of some interest to rcderivc the limiting value of du/dv by using other generic relationships between u 
and V. As before, to avoid 0/0 forms, we work with Sectors III & IV. In particular, in Sector III, we have 

- = coth (45) 
V Am 

By differentiating this equation w.r.t. v, we obtain 

1 du u 1 1 dt 



V dv 4msinh^(t/4m) dv 

By recalling that sinh(i/4TO) = v/fi, we rewrite the above Eq. as 

du u _ —f1 1 dt 

dv V Am V dv 

Now, from Eqs. (10) and (38), note that 

dt dt dr AmE 



dv dr dv Vy/E'^ - 1 + 2m/r - uE 

And the limiting value of 



(46) 



(47) 



(48) 
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dv u ~ V uh + uh uh 
And since /i(2to) = 0, we find from Eq. (47) that 

du u 



dt Am Am 2m , , 

— r ^ 2m (49) 



= 0; r = 2m (50) 

dv V 



Or, 

du u Uh 



r = 2m {Sector III + IV) (51) 



dv V —Uh 

Similarly, for the sake of overall consistency, in Sectors, I & III, we must have du/dv = +1 at r = 2m. 



V. A DIFFERENT CONSIDERATION 



Actually we could have obtained the above derived unique result in a relatively simpler manner by differentiating 
the Global Eq. 

^2 - ^2 = (r/2TO - l)e'^/2™ (52) 

w.r.t. V : 

2u^-2v= ^ (53) 

dv Am'' dv 

First let us note from Eq.(37) that in Sectors II & IV, the hmiting value of 

dr 4m(l — 2m/r) 2m(l — 2m/r) 
dv V — u Vh 

Then, by using this above Eq. in (53), we find 



0; r ^ 2m (54) 



du „ 
Uh = vh; r ^ 2m (55) 



so that 

du Vh 



dv —Vh 



-1; r^2m (56) 



VI. CONCLUSIONS 



The Kruskal coordinates were found way back in 1960, and in the present paper, we have worked out some aspects 
of the kinematics of a test particle following a radial Kruskal geodesic. To attain this we used, for the first time, the 
precise value of uh and Vh as a function of the initial conditions of the problem ri, m or E. It is clearly found that 
u^j = v^ is non-zero in general. 

We then proceeded to obtain the expression for the Kruskal derivative in terms of m, E and r. We found that the 
Kruskal derivative is regular at the EH unlike the Scharzschild deivative(s) where dt/dr = — oo at the EH. 

In particular du/dv — +1 at the EH if we consider the "other universe" whose existence is suggested by the full 
Kruskal manifold, and which is a time reversed version of "our universe". But, if we move to the "our universe", the 
expected value of du/dv = 1 at the EH. 

The regular nature of the Kruskal derivative is in keeping with the notion that Kruskal coorinates are free of 
singularities at the EH. 

In a subsequent paper, we shall find out other important features of the Kruskal dynamics vis-a-vis the well known 
Schwarzschild dynamics. 
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